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J A  random  vibration  predictive  model  for  cable  strumming  response  to 
ocean  currents  is  presented.  The  model  is  justified  on  the  basis  of  an 
extensive  examination  of  four  different  field  experiments,  conducted  for 
the  purpose  of  studying  vortex  induced  cable  vibration.  The  examination 
of  these  data  reveals  that  real  world  cable  strumming  is  predominantly 
a  random  process  that  cannot  be  predicted  on  the  basis  of  extending  lab¬ 
oratory  observations  of  essentially  deterministic  phenomena  such  as 
lockin.  Representative  data  are  presented  to  support  these  conclusions. 

The  proposed  predictive  model  is  demonstrated  by  two  examples:  the 
response  of  a  uniform  cable  in  a  non-uniform  flow  and  the  response  of  a 
cable  non-uniform  in  diameter  to  a  uniform  flow.  The  results  of  the 
second  example  are  favorably  compared  to  measurements  made  in  a  simple 
field  test.^  Z 

INTRODUCTION 

The  prediction  of  cable  vibration  response  to  incident  currents, 
cable  strumming,  is  an  as  yet  unresolved  problem  in  cable  systems 
design. 

Observations  of  real  world  cable  systems  reveal  the  phenomena  to  be 
both  complex  and  quite  different  from  published  laboratory  results.  The 
reasons  for  this  are 'that  field  conditions  rarely  possess  the  temporal 
and  spatial  uniformities  of  the  currents  used  in  laboratory  tests,  that 
long  field  cables  have  much  higher  modal  densities  at  the  frequencies  of 
interest,  and  that  field  cables  are  rarely  uniform  in  properties  over 
their  length,  as  contrasted  to  the  short  uniform  sections,  commonly 
chosen  for  laboratory  work. 

In  particular,  the  exhaustively  studied  phenomenon  known  as  reson¬ 
ant  lockin,  or  synchronization  of  the  wake  with  a  natural  frequency,  is 
a  relatively  uncommon  aspect  of  the  response  of  field  deployed  cables. 
Instead  of  the  deterministic,  periodic  response,  characteristic  of 
lockin,  the  motion  time  history  of  a  typical  offshore  cable  system  is 
usually  best  described  as  a  random  process.  The  periodic  response  of  a 
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single  mode  is  quite  rare;  more  often  the  response  is  a  superposition  of 
from  a  few  to  hundreds  of  modes  responding  simultaneously. 

The  purpose  of  this  paper  is  to  propose  a  simple  stochastic  model 
for  the  prediction  of  the  non-lockin  cable  vibration  response  so  fre¬ 
quently  observed  in  the  field.  The  model  was  formulated  after  an  in 
depth  examination  of  the  data  collected  on  four  separate  and  quite  dif¬ 
ferent  field  experiments.  The  individual  experiments  and  principal  con¬ 
clusions  are  reviewed  to  establish  a  physical  basis  for  the  model.  The 
model  is  presented  and  its  principal  features  are  demonstrated  in  two 
example  calculations.  In  the  first  example  the  response  of  a  uniform 
cable  to  a  spatially  non-uniform  flow  is  predicted.  The  second  example 
is  the  response  prediction  for  a  cable  non-uniform  in  diameter  over  its 
length  exposed  to  a  uniform  flow. 

A  field  experiment  is  described  in  which  the  measured  response  of  a 
cable  non-uniform  in  diameter  is  compared  to  the  response  predicted  in 
the  example. 

EXPERIMENTAL  EVIDENCE 

Laboratory  Tests 

The  literature  concerned  with  fixed  and  moving  cylinders  is  exten¬ 
sive.  The  principal  results  are  presented  in  a  recent  review  paper  by 
King  [3].  Laboratory  tests  of  cables  are  not  so  numerous  and  have  been 
generally  restricted  to  Investigation  of  cross  flow  lockin  phenomena  on 
relatively  short  lengths  of  cable  (0.5  to  5  meters).  The  flows  have 
been  typically  low  in  turbulence  and  spatially  uniform.  Only  the  first 
few  natural  response  modes  have  been  studied  and  then  only  one  at  a  time. 
Some  valuable  conclusions  are  that  at  lockin  the  response  is  self  limit¬ 
ing  with  the  antinode  responses  approaching  one  diameter  single  ampli¬ 
tude.  Shedding  at  the  antinodes  is  similar  to  locked-in  moving  cylin¬ 
ders  and  shedding  at  the  nodes  is  like  that  of  non-moving  cylinders  [8]. 
References  [1]  and  [6]  are  additional  sources  of  experimental  cable 
strunming  results. 

The  laboratory  experiments  have  made  valuable  contributions  to  our 
understanding  of  the  local  interaction  between  the  motion  of  the  cylin¬ 
der  and  the  wake,  but  cannot  simulate  the  high  modal  densities  and 
spatial  and  temporal  variations  in  the  flow.  The  response  of  actual 
field  deployed  cables  must  be  examined  to  reveal  these  effects. 


Field  Tests 

Four  field  experiments  are  reviewed  below.  The  raw  data  from  three 
of  them  were  reduced  and  analyzed  by  the  authors  in  the  past  18  months. 
The  previously  reduced  data  from  the  fourth  experiment  known  as  "fish- 
bite"  were  also  evaluated  and  reported  here. 

Bermuda  Testspan:  From  December  1973  to  February  of  1974  the 
United  States  Navy  conducted  a  cable  strumming  experiment  on  Plantagent 
Bank  off  Argus  Island,  Bermuda.  A  256  meter  long,  .016  m  diameter 
electromechanical  cable  was  stretched  horizontally  at  a  depth  of  27.5 
meters.  The  cable  without  strumming  suppression  devices  was  mechan¬ 
ically  very  complex.  Numerous  instrument  pods,  lead  weights,  and  floats 
were  distributed  over  its  length.  Two  Savonlous  rotor  current  meters 
suspended  near  mid-span  revealed  extreme  variation  in  both  current  speed 
and  direction.  Flow  reversals  occurred  with  an  average  period  of  eight 
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seconds,  suggesting  that  long  wavelength  surface  waves  were  influencing 
the  flow.  The  flow  speed  was  generally  less  than  0.5  m/s.  Instrumen¬ 
tation  failures,  signal  to  noise  problems,  and  inadequate  accelerometer 
calibration  data  made  the  analysis  of  the  five  year  old  tapes  quite 
difficult.  Nonetheless,  the  following  conclusions  were  reached. 

The  cable  strumming  response  was  typical  of  a  broad  band  random  process. 
Mode  participation  varied  with  flow  speed  from  the  10th  to  the  150th 
mode.  Figure  1  shows  a  typical  acceleration  response  spectrum  at  a 
point  on  the  cable.  This  spectrum  represents  the  sum  of  the  two  spec¬ 
tral  components  from  a  biaxial  accelerometer.  Due  to  cable  rotation  it 
was  not  feasible  to  separate  cross  flow  from  in  line  cable  response. 

This  particular  spectrum  was  averaged  from  15,  80-second  records.  It 
reveals  significant  (within  lOdB  of  the  peak)  response  over  the  range  of 
2  to  4  Hz.  More  than  10  modes  are  simultaneously  contributing  to  the 
response.  The  conclusions  drawn  from  examination  of  many  hours  of  data 
were  that  a)  lockln  was  rare,  b)  the  response  was  generally  broadband 
random,  and  c)  over  many  hours  the  root  mean  square  (rms)  response  was 
self-limiting  and  constant  (to  within  2  dB) .  Though  suspect  because  of 
poor  calibration  data,  the  rms  response  estimate  was  less  than  a  cable 
diameter. 

It  is  the  authors’  conclusion  that  the  combination  of  rapid  varia¬ 
tions  in  flow  velocity,  non-uniform  cable  properties  (especially  dia¬ 
meter)  and  high  modal  density  are  responsible  for  the  broadband,  multi- 
moded  response  observed  in  this  test.  Additional  description  of  this 
experiment  may  be  found  in  reference  [21. 

Seacon  II:  The  experiment  was  conducted  by  the  Naval  Civil  Engin¬ 
eering  Laboratory  in  1975  [4,9].  It  consisted  of  a  horizontal  delta 
array  buoyed  and  anchored  at  three  corners,  152  meters  below  the  surface 
in  880  meters  of  water  in  the  Santa  Monica  Basin.  One  of  the  horizontal 
arms  of  the  delta  was  instrumented  with  accelerometers  to  measure  the 


cable  strumming  response.  The  unfalred  cable,  sensors,  and  electronics 
were  essentially  the  same  as  that  used  in  the  Bermuda  experiment.  No 
lead  weights  or  buoyancy  bags  were  used. 

In  this  experiment  the  current  speeds  were  on  the  order  of  0.1  m/s, 
were  tidally  generated,  and  were  slowly  varying.  The  current  meters 
were  more  than  200  meters  from  the  test  cable  and  the  current  and  res¬ 
ponse  data  were  not  well  synchronized.  Data  reduction  and  analysis 
difficulties  encountered  in  evaluating  the  four  year  old  raw  data  tapes 
were  similar  to  those  encountered  with  the  Bermuda  data.  Due  to  the 
lower  current  speeds  the  participating  natural  modes  ranged  from  the 
2nd  to  the  15th  with  from  two  to  six  modes  responding  simultaneously. 

The  vector  averaged  (cross  plus  in  line)  rms  response  was  self  limiting, 
lockln  was  rare,  and  multi-moded  random  vibration  behavior  predominated. 
Even  though  the  current  was  relatively  uniform  the  response  typically 
spanned  a  considerable  range  in  frequency.  It  was  concluded  that  this 
was  a  consequence  of  a  large  range  of  cable  and  instrument  housing  diameters 
(from  .016  to  0.15  m).  Fot  additional  Information  consult  reference  [2]. 

Fishbite:  Fishbite  is  the  name  given  to  a  long  cable  experiment 
conducted  in  1976  by  Softley,  Dilley  and  Rogers  [10]. 


15  HZ 


FIGURE  2.  FISHBITE:  RMS  DISPLACEMENT 
SPECTRUM 

A  .012  m  diameter  wire  rope  500  meters  long  was  hung  over  the  side 
of  a  ship  anchored  in  1960  meters  of  water  at  The  Tongue  of  the  Ocean, 
77s  52 'W  and  25*  10’N.  The  tidal  flow  was  neither  spatially  nor  tem¬ 
porally  uniform,  varying  from  0.1  to  0.4  meters  per  second.  The  cable 
was  of  uniform  diameter  with  a  current  meter  attached  at  the  halfway 
point.  No  other  lumped  masses  were  attached  to  the  cable.  The  response 
was  measured  at  the  top  end  only.  The  modal  spacing  was  .025  Hz. 

The  repsonse  typically  Included  more  than  a  hundred  modes  between 
8  and  12  Hz,  centered  on  10  Hz.  A  root  mean  square  displacement  spec¬ 
trum  (from  Softley  et  al..  Figure  28)  is  shown  in  Figure  2.  The  rms 
response  at  the  upper  end  was  limited  to  less  than  a  diameter.  The  band 
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width  o£  the  response  can  only  be  attributed  to  the  spatial  and  temporal 
variation  of  the  current.  The  pertinent  conclusions  are  that  the  re¬ 
sponse  was  rms  limited  and  broadband  random.  Lockln  was  not  observed. 

Castine,  Maine:  Field  tests  were  conducted  near  Caatine,  Maine, 
during  the  summers  of  1975  and  1976  by  Vandiver  and  Mazel  [5,  7,  11]. 
Figure  3  is  a  schematic  illustration  of  the  test  site,  which  was  located 
on  a  sandbar  that  was  flooded  at  high  tide  and  exposed  at  low  tide. 

Tests  were  conducted  on  the  rising  tide.  The  flow  was  normal  to  the 
23.3  m  long  cable  and  varied  from  a  maximum  of  0.7  m/s,  shortly  after 
submergence,  to  0  m/s  at  high  tide  over  a  period  of  approximately  2-1/2 
hours.  The  flow  was  spatially  uniform  to  within  5Z  over  the  length  of 
the  cable.  The  mean  flow  speed  variations  could  be  considered  as  quasi¬ 
static  when  viewed  on  a  time  scale  of  many  periods  of  cable  vibration. 

The  data  were  taken  in  calm  weather  and  surface  waves  were  not  a  problem. 


High  Tide  Level  (3m  deep) 


Cable 

Tensiometer  Accelerometer  \  winch 


FIGURE  3.  CASTINE  EXPERIMENT  SET-UP 


The  instrumentation  consisted  of  a  tensiometer,  an  electro-magnetic 
current  meter,  very  small  moveable  accelerometers,  and  a  direct  dis¬ 
placement  measuring  device,  called  the  "fish".  The  fish  and  accelero¬ 
meters  could  measure  cross  flow  vibration  without  loading  the  cable. 

Many  different  types  of  cables  with  and  without  strumming  suppression 
devices  were  tested  during  two  months  of  on-site  testing.  The  cables 
tested  were  of  various  diameters  from  .006  to  .016  m. 

The  cable  strumming  was  all  in  the  subcritical  Reynolds  number 
range,  and  the  vibration  response  typically  included  one  or  more  of  the 
first  fifteen  natural  modes.  Vibration  in  line  with  the  flow  was  negli- 
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gible  when  compared  to  the  cross  flow  vibration. 

The  experiments  were  conceived  as  a  stepping  stone  from  idealized 
laboratory  conditions  to  the  almost  hopelessly  complex  conditions  typical 
of  the  Bermuda  test  and  Beacon  IX.  Deterministic  lockin  response  and 
non-deterministlc  non-lockin  were  observed  at  Castine. 

Single  frequency  sinusoidal  lockin  response  occurred  at  discrete 
steps  in  flow  speed  for  all  cables  tested.  Responding  modes  varied  from 
the  first  to  as  high  as  the  15th,  depending  on  the  particular  cable 
properties,  tension,  and  flow  velocity.  The  Reynolds  number  range  for 
observed  cross  flow  vibration  was  from  300  to  8000.  Antinode  response 
was  limited  to  approximately  one  diameter  for  all  cables  tested.  In 
general  lockin  behavior  was  exhibited  over  a  smaller  range  of  reduced 
velocity  (U_  ■  U/f  D)  than  observed  in  the  laboratory  with  vibrating 

rigid  cylinders.  U  ‘’is  the  flow  velocity,  f  is  the  frequency  of  the 
largest  peak  in  the  response  spectrum,  and  D  ^  is  the  cable  diameter. 
The  extremes  of  the  lockin  range  are  evidently  reduced  by  real  world 
spatial  and  temporal  variations  in  the  flow,  not  revealed  by  a  single 
point  measurement  of  mean  flow  speed.  As  the  mean  flow  speed  was  ob¬ 
served  to  slowly  move  away  from  a  value  that  insured  lockin,  the  response 
became  progressively  more  random.  This  suggests  that  the  portion  of  the 
cable  near  the  antinodes  experiencing  lockin  was  decreasing. 

Even  at  flow  speeds  which  completely  prevented  wake  synchronization 
with  a  natural  frequency  of  the  cable,  the  observed  response  was  quite 
substantial.  The  observed  response  at  Castine  under  such  conditions 
was  best  described  as  a  narrow  band  random  process.  The  term  "non- 
lockin"  was  coined  to  describe  this  and  all  other  random  vibration  cable 
response.  Non-lockin  response  occurred  about  half  of  the  time  at  Cas¬ 
tine,  and  as  mentioned  before  dominated  the  responses  observed  in  the 
Bermuda,  Seacon  II  and  fishblte  experiments. 

As  a  general  classification  non-lockin  occurs  when  variations  in 
the  current  or  non-uniformities  in  the  cable  prevent  widespread  wake 
synchronization  with  a  single  natural  frequency  of  the  cable. 

It  is  conjectured  that  under  such  conditions  the  correlation  length 
of  shed  vortices  is  on  the  order  of  only  a  few  diameters  and  that  the 
motion  of  the  cable  at  any  point  is  poorly  correlated  to  the  wake 
behind  it.  Under  such  conditions  the  lift  force  spectrum  at  a  point  on 
the  cable  can  be  characterized  as  a  band  limited  random  process  whose 
center  frequency  is  adequately  characterized  by  an  empirically  deter¬ 
mined  constant  such  as  a  reduced  velocity  or  a  Strouhal  number.  In  many 
respects  this  proposed  model  of  a  lift  force  spectrum  is  similar  to  that 
measured  on  stationary  cylinders. 

Spectra  of  non-lockin  displacement  data  from  Castine  typically 
revealed  one  to  four  closely  grouped  natural  frequencies  of  the  cable, 
and  occasionally  the  shedding  frequency  which  would  be  expected  for  a 
stationary  cylinder  in  that  flow.  The  narrow  banded  nature  of  the  non- 
lockin  response  observed  at  Castine  is  consistent  with  the  test  condi¬ 
tions:  uniform  diameter  cables  in  a  quite  uniform  flow.  The  non-lockin 
rms  cross  flow  response  was  approximately  one  quarter  of  a  diameter  for 
all  cables  tested  and  was  not  measurably  Reynolds  number  dependent. 

This  response  level  is  considerably  less  than  lockin  response  amplitudes. 
However,  because  the  displacements  are  a  random  process,  individual 
peaks  often  exceeded  two  diameters.  The  displacements  could  be  ade¬ 
quately  described  by  a  Gaussian  random  process.  When  measured  in  terms 
of  diameters,  the  rms  response  was  relatively  insensitive  to  cable 
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material,  tension,  flow  velocity,  and  the  number  of  participating  nodes. 
Furthermore,  the  reduced  velocity  remained  nearly  constant  for  each  type 
of  cable  tested  for  all  flow  conditions.  The  observed  mean  values  for 
reduced  velocities  varied  among  the  cables  from  4.5  to  6.7.  Additional 
data  from  Castlne  will  be  presented  later  in  this  paper. 

At  the  end  of  this  extensive  review  of  field  data  it  was  concluded 
that  a  predictive  model  for  non-lockin  cable  vibration  response  was 
possible  and  necessary.  Such  a  model  would  be  useful  for  engineering 
design  purposes  if  it  could  provide  an  estimate  of  the  response  spectrum 
of  a  deployed  cable.  From  the  available  experimental  data  it  was  felt 
that  such  a  predictive  model  could  be  based  upon  the  concept  of  a  self- 
limiting  response  and  a  band  limited  local  lift  force  spectrum  centered 
at  a  frequency  established  by  the  local  flow  velocity  and  cable  diameter. 

A  STOCHASTIC  MODEL 

The  review  of  field  data  briefly  described  in  the  previous  section 
indicates  that  deterministic  lockin  models  are  rarely  appropriate  in  the 
prediction  of  real  world  cable  strumming  response.  Hon- lockin  random 
vibration  generally  dominates.  In  this  section  a  simple  stochastic 
model  is  proposed  for  the  prediction  of  random  vibration  cable  response 
to  current  flows. 


In  this  analysis  the  linearizing  assumption  is  made  that  the  cross 
and  in  line  cable  excitations  and  responses  can  be  decoupled.  Since  the 
cross  flow  response  is  typically  an  order  of  magnitude  greater  than  the 
in  line  response  it  is  the  only  response  considered  here. 

The  quantities  of  interest  in  a  stochastic  framework  are  mean 
squared  responses  and  power  flows.  In  cable  strumming  the  spatially 
distributed  lift  forces  arising  from  the  local  vortex  shedding  along  the 
cable,  constitute  the  power  input  into  the  cable  system.  Qualitatively 
and  quantitatively  very  little  is  known  about  the  mechanism  that  governs 
the  input  power.  It  appears  to  be  a  non-linear  feedback  based  on  mean 
square  amplitude.  For  the  purpose  of  this  model  it  is  assumed  that  the 
feedback  mechanism  is  adequately  modelled  by  a  uniform  spatially  distri¬ 
buted  feedback  level  based  on  the  cable's  spatially  averaged  mean  square 
displacement  response.  The  cable  length  is  finite  and  the 
principal  power  outputs  are  due  to  internal  cable  friction  and  fluid 
damping.  As  the  system  is  spatially  distributed,  power  is  shared 
throughout  the  cable  and  power  output  at  any  point  is  drawn  from  the 
entire  cable  system. 

With  these  simplifications  a  form  of  the  spatially  distributed 
excitation  may  be  proposed.  f(x,t)  is  the  force  per  unit  length  at  a 
location  x  on  the  cable  at  time  t  and  is  given  by 

f(x,t)  -  .5pwD(x)U2(x)CL(x,t)  (1) 

Py:  water  density 
D(x):  cable  diameter  at  x 

U(x):  normal  component  of  the  freestream  velocity  at  x 
CL(x,t):  local  momentary  lift  coefficient 

The  local  lift  coefficient  CL(x,t)  is  a  random  variable  which  is 


assumed  to  be  zero  mean  Gaussian,  time  statlonary-ergodlc,  and  locally 
independent  of  cable  response.  The  first  two  assumptions  are  straight¬ 
forward  and  reasonable.  The  third  arises  from  the  previous  statement 
that  a  uniform  feedback  level  is  obtained  from  a  spatially  averaged  mean 
square  response.  As  will  be  demonstrated  this  assumption  is  acceptable 
for  cables  of  uniform  diameter  in  non-uniform  flows,  but  presents  obsta¬ 
cles  for  cables  with  non-uniform  diameters. 

Given  the  above  assumptions  we  can  completely  characterize  CL<x,t) 
by  either  its  correlation  function  or  its  power  spectrum. 

CORRELATION:  Rc<x,y,T)  -  EfC^x.t)  (^‘(y.t+x)]  (2) 

/*»  i2wfT 

dT  Rc(x,y,T)e  (3) 

From  experimental  results  we  know  that  the  spatial  correlation 
length  for  non-lockln  is  on  the  order  of  a  few  cable  diameters.  If  the 
assumption  is  made  that  the  shortest  excited  spatial  wavelength  is  much 
larger  than  this  correlation  length,  then  the  lift  coefficient  reduces 
to: 


S  (x,y,f)  *  L6(x  -  y)S  (x,f)  (A) 

C  v 

L:  cable  length 

6(x-y):  impulse  function  at  x  *  y 

Because  the  units  of  6(x-y)  are  length  1 ,  L,  the  length  of  the  cable, 
is  introduced  in  the  above  equation  to  maintain  consistent  units.  A 
visual  image  of  this  spatially  incoherent  process  is  that  of  "rain  on 
the  roof”. 

A  box  car  model  for  the  spectrum  of  the  lift  coefficient  is  pro¬ 
posed  and  illustrated  in  Figure  A. 


f  (x)  -  S  !?{*?■  ,  local  center  frequency  of  (5) 

8  vortex  shedding 

W(x)  *  2P  f  (x)  ,  bandwidth  of  shedding  (6) 

frequencies 
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c 


2WH 


Sc(*,f)d£ 


the  spatially  uniform  (7) 

mean  square  lift  coefficient 


This  model  for  the  lift  coefficient  spectrum  is  dependent  upon 
three  dimensionless  parameters: 

1)  S£:  the  cable's  non-lockin  Strouhal  number.  It  determines 

the  center  of  the  spectrum  and  is  inferred  from  ex¬ 
perimental  observation  of  the  reduced  velocity. 

2)  Pt :  the  excitation  bandwidth,  as  a  percentage  of  the 

S£  dependent  shedding  frequency. 

3)  C:  mean  square  lift  coefficient  appropriate  to  the 

cable,  flow  conditions,  and  feedback  level. 


C  is  spatially  distributed  and  uniform,  and  incorporates  all  cable- 
fluid  feedback  mechanisms.  Its  value  will  be  set  such  that  the  spttially 
averaged  mean  square  response  is  at  its  appropriate  value.  A  more  exact 
form  for  Sc(x,f)  would  include  tails  at  higher  and  lower  frequencies. 

The  boxcar  model  is  felt  to  be  appropriate  for  the  simple  response  pre¬ 
diction  model  developed  here.  Refinements  may  be  warranted  in  the 
future. 

This  simple  model  for  the  mean  square  lift  coefficient  may  be  used 
with  equation  (1)  to  derive  an  expression  for  the  distributed  lift  force 
spectrum. 


Sp(x,y,f)  - 


/ 


dTE[f(x,y)f(x,t+T)Je 


12irfT 


(8) 


(.5pw)2  (D(x)U(x))3  CL6(x-y)B(f,x) 

4P  S 
t  t 


(9) 


B(f,x)  -  1 

-  0 


(1-Vst  5$ s  l£l  <  (lwt>s,  5oo 

elsewhere 


(10) 


The  cable  model  must  be  obtained  next.  The  key  assumption  is  that 
it  is  modelled  as  a  linear  time  invariant  system.  A  modal  decomposition 
approach  is  used  here.  This  requires  knowledge  of  the  individual  natu¬ 
ral  frequencies  f  ,  mode  shapes  ^ (x) ,  and  modal  damping  ratios  . 

The  modal  transfer  functions  Hj (f )  for  a  cable  with  fixed  ends,  and 

virtual  mass  per  unit  length  m(x)  are  completely  specified  by: 


(ID 


y  dx  ♦j(x)4>k(x)  m(x) 
0 


(12) 


9 


(13) 


Hj  (x) 


+  2U, 


The  nodal  force  spectrum  la  obtained  from  the  node  shapes  and  dis¬ 
tributed  force  spectrun  as  shown  In  the  next  equation. 


L  L 

J*  dx  J  dy  ^(x)  $fc(y)  SF(x,y,f) 


In  this  analysis  the  contributions  to  the  response  resulting  from  the 
modal  overlap  cross  terns  (j  j*  k)  are  assumed  snail  compared  to  the 
terms  for  which  j  -  k.  Damping  generated  modal  coupling  Is  also  neglec¬ 
ted.  With  these  assumptions  the  response  spectrum  may  be  obtained  as 
follows : 

SD(x,y,f)  “^(x)  ♦j(y)lHj(f)|2  SF^(f)  (15) 

where  S„  (f)  Is  reduced  from  equation  (14) ,  and  takes  the  form: 

J  rL  rL 

SF  (f)  -  J  dx  J  dy  4»j(x)()>j(y)SF(x,y,f)  (16) 

J  0  0 

L 

(.5pw)2  CL  f  dx  fj2  (x)^D(x)U(x)J3B(f,x)  (17) 

4PtSt  * 

The  total  displacement  response  spectrum  is  therefore  the  sum  of  the 
Individual  modal  response  spectra,  summed  over  all  modes  included  in 
the  excitation  bandwidth  defined  by  the  function  B(f,x)  from  equation 
(10). 

It  remains  to  establish  values  for  the  three  parameters  C,  and 

P  .  C  results  from  an  empirically  Imposed  limit  on  the  spatially  aver- 
C  2 

aged  mean  square  response  A  ,  expressed  in  diameters  squared,  where 


A2  -  f  /L dx  4-  /df  S  (x,x,f) 
L  0  D 2  D 


In  other  words,  C  is  adjusted  so  that  the  empirically  observed  limits  on 
A,  the  rms  response,  are  satisfied.  Given  experimentally  determined 
values  for  A,  and  St,  the  response  is  now  completely  characterized. 

Modal  Parameter  Values 

The  data  from  six  cables  tested  at  Castine  were  used  to  develop  a 
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preliminary  data  base  for  A,  Pfc  and  Sfc.  The  non- lock in  cross  flow  dis¬ 
placement  time  histories  were  divided  into  five  second  records  from 
which  individual  raw  spectra  were  computed.  The  root  mean  square  dis¬ 
placement  obtained  from  each  spectrum  when  divided  by  the  appropriate 
cable  diamete^  provided  an  estimate  of  A.  S  was  estimated  as  the  reci¬ 
procal  of  the  observed  reduced  velocity 


where  f  :  frequency  of  the  highest 
P 

dividing  the  5  dB  down  bandwidth  of 
center  frequency  of  the  band. 


spectral  peak.  P  was  estimated  by 
the  response  spectrum  by  twice  the 


P  *  Af5dB  (20) 

2f 

P 


Five  dB  down  was  an  arbitrary  definition  of  significant  response  level. 
For  some  cases  the  response  spectrum  contained  only  one  sharp  modal  re¬ 
sponse  peak.  In  such  cases  only  an  upper  bound  estimate  on  P  could  be 
estimated  by  dividing  the  computed  modal  separation  frequency  by  twice 
the  frequency  of  the  observed  peak.  Table  1  briefly  describes  the  six 
cables.  Material,  diameter,  virtual  mass  per  unit  length  and  the  ten¬ 
sion  range  are  given.  The  flow  speed  varied  with  the  tidal  cycle.  The 
virtual  mass  per  unit  length  m,  was  estimated  by  adding  the  displaced 
mass  per  unit  length  of  water  (the  added  mass)  to  the  wet  (saturated) 
mass  per  unit  length  of  the  cable.  Table  II  presents  the  results  of  the 
parameter  estimation  for  A,  S  and  Pt.  The  nunfter  of  records  used  and 
the  standard  deviation  for  each  value  are  given.  Upper  bound  estimates 
are  denoted  by  the  note  U.B.  in  the  standard  deviation  column. 

EXAMPLE  RESPONSE  PREDICTIONS 


A  Uniform  Cable  in  a  Non-Uniform  Flow 

A  long  cable  of  length  L  and  diameter  D  with  fixed  ends  and  con¬ 
stant  tension  is  exposed  to  a  normal  flow  field,  which  is  constant  but 


different  for  the  two  portions  of  the  cable. 

U(x)  -  Ux  0  <  x  <  Lx  (21) 

-  U2  Lx  <  x  <  L  (22) 

The  natural  frequencies  and  mode  shapes  are  specified  by: 
fj  -  JCp/2L  (23) 

4>j  (x)  *  /2  sin  (24) 


where  Cp:  phase  velocity  of  transverse  waves. 

By  assuming  that  the  cable  is  long  compared  to  the  longest  excited 
spatial  wavelength,  a  piecewise  approximation  to  the  integral  contained 
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Table  I.  CASTINE  TEST  CABLES 


Description 

Diameter 

(m) 

Virtual  Mass 
per  unit  length 
(kg/m) 

Range  of 
Tension 
(nt) 

Sampson,  "Blue  Streak",  12  strand 
single  braid,  polyester  and  poly¬ 
propylene 

.010 

.177 

310-1025 

Philadelphia  Resins  "Philly- 
stran",  7  strand  Kevlar  core 
rope  with  a  polyurethane  jacket 

.0123 

.233 

490-2000 

Wire  rope,  3x9  torque  balanced 
galvanized  plow  steel  with  a 
polyethylene  jacket 

.0071 

.148 

270-2600 

Wall  Rope  Works  "Uniline" 
polyester  cable,  woven  jacket 

.0131 

.290 

1150-1500 

Philadelphia  Resins,  PS  29  EM  1 
Four  conductors  inside  single 
braided  Kevlar,  protected  by  a 
woven  Dacron  jacket 

.0071 

.110 

850-1110 

Cortland  Line  Co.,  parallel 

.0061 

.081 

1700-2200 

fiber  Kevlar  with  a  braided  jacket 

************************* 


Table  II.  PARAMETER  ESTIMATES  (MEAN 
AND  STANDARD  DEVIATION) 


Cable 

Number  of 
Samples 

A 

°A 

St 

°S 

Pt 

°P 

Blue  Streak 

140 

.26 

.04 

.15 

.020 

.08 

.03 

Phillystran 

43 

.20 

.04 

.22 

.015 

.10 

.03 

Wire  Rope 

6 

.28 

.07 

current  meter  failed 

Uniline 

12 

.20 

.02 

.19 

.004 

.10 

U.B.* 

PS  29 

40 

.28 

.05 

.16 

.007 

.12 

U.B.* 

Cortland 

43 

.29 

.04 

.17 

.010 

.13 

U.B.* 

^Denotes  upper  bound 


in  equation  (17)  may  be  obtained. 


/Ldx  <)>j2(x)(D(x)U(xi3B(f,x)  *  D3[U13L1B1(f)  +  U23L2B2<f)]  (25) 

where  B^f)  -  1  (l-Pt)St  ^  <  |f |  <  (l+Pt)St  ^ 

■  0  elsewhere  (26) 

»2<f>  ■  1  (1-pt>st  If  I  <  <l+Pt>St  T 

*  0  elsewhere  (27) 


This  leads  directly  to  an  expression  for  the  modal  exciting  force  Sp  (f) 
and  hence  to  the  response  spectrum.  j 


2  3 

p  CLD 

SD(x,y,f)  g  16P~S~  ^  ♦J<*)#J(y)|"J(f)r 


[U13L1B1(f)+U3L2B2(f)]  (28) 


The  relative  impact  of  the  flows  IL  and  U2  on  the  response  may  now  be 
evaluated.  Let  P  (x)  and  P2(x)  be  the  local  mean  square  responses  at 
a  point  x  due  to  the  two  flow  fields.  P  (x)  and  P2(x)  may  be  obtained 
by  requiring  that  x  *  y  and  integrating  the  expression  for  Sp(x,x,f) 
over  the  frequency  bands  defined  by  B^(f)  and  B2(f)  respectively.  The 
result  for  Pj(x)  is  shown  in  equation  (29) .  The  result  for  P^(x)  may  be 
obtained  by  replacing  the  subscript  1  by  a  2  in  equation  (29)7 


Px(x) 


/ 


df 


2  3 

p  CLD 
w 


16 


P  S 
t  t 


s 

i 


<|»j2(x)|Hj(f)|2  U13L1B1(f) 


(29) 


P^(x)  and  P2(x)  depend  on  the  mode  shapes  (x) .  To  eliminate  this  x 
dependence,  spatially  averaged  values  may  be  computed,  denoting  the 
result  as  P^  and  P2-  For  example,  P^  is  given  by 


P 


1 


00 

/  df 


2_  1 
p  CLD 
w _ 

16  PtSt 


£^|Hj(f) 


|2  U1\1B1(f) 


(30) 


The  relative  spatially  averaged  mean  square  response  is  given  by 
£j  |HJ(t)|2B1(« 

'2  U23l2  f*  l»1<f)|2B2<t) 

J  *00  J 


The  integration  over  frequency  can  be  simplified  by  the  following 
approximation: 


1) 


The  modal  damping  ratios 


were  assumed  small  (less  than  152) . 


A 
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mode  whose  natural  frequency  was  inside  the  B. (f)  or  B.(f)  bands 
was  assumed  to  have  its  modal  bandwidth  completely  inside  B. (f)  or 
B2(f).  1 

2)  The  spatially  averaged  mean  square  response  of  any  mode  inside 

B^(f)  or  B2(f)  was  assumed  equal  to  any  other  mode  Inside  B1(f)  or 
B2(f)  respectively. 

Approximation  1  allows  a  simple  expression  for  the  response  of  a 
single  mode  j  to  a  white  noise  force  spectrum  of  unit  amplitude. 


H.  (f)  1 2df  -  - - - T  ‘  (32) 

3  n  2Cj(2irfj)J 

Approximation  2  allows  the  response  of  each  mode  included  in  the 
band  to  be  represented  by  the  response  of  an  equivalent  mode  which  has 
a  natural  frequency  equal  to  the  center  frequency  of  the  band  and  a 
damping  ratio  typical  of  all  modes  in  the  band.  Therefore,  within  each 
band  the  following  relations  apply. 


u. 

Vf>! 

*  5.  -  constant,  f. 

^1  3 

s  if 

t  D 

(33) 

B2(f)t 

*  5,  -  constant,  f. 

h  3 

\  ' 

U2 

S  — 
t  D 

(34) 

The  number  of  modes  in  a  band  is  the  product  of  the 
modal  density.  For  a  uniform  taut  string  the  modal 
the  reciprocal  of  the  fundamental  natural  frequency 
modes  per  Hz.  Therefore  the  number  of  modes  inside 
given  by 

2P  S  (35) 

D 

2P  S  ^2  (36) 

D 


bandwidth  and  the 
density  is  simply 
f.  and  has  units  of 
B^(f)  and  B2(f)  are 


The  total  response  of  all  modes  in  the  band  is  given  by  the  product  of 
the  number  of  equally  responding  modes  and  the  response  of  the  typical 
mode. 

For  example,  for  the  band  B^(f) 


£ 
i  . 


/ 


df |H  (f)| 


Bx(f) 


2L 

C 

P 


2PtSt  D 


_ 1 _ 

25.  M2  (2irf  )3 
J1  J1 


A  similar  result  is  obtained  for  the  band  B2(f)  by  replacing  the  sub* 
scripts  1  with  a  2.  These  results  may  now  be  substituted  in  equation 
(31)  to  obtain  Pj/P2*  'Hie  result  depends  on  the  choice  of  damping 
model  assumed . 

Two  plausible  models  are: 
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Constant  damping  ratio 

£.  "  S4 

J1  J2 

Frequency  dependent  damping  ratio 


fj 

J1  J1 


f j 

J2  J2 


(38) 

i 


(39) 


A  frequency  dependent  damping  ratio  la  equivalent  to  specifying  a  con¬ 
stant  damping  coefficient  per  unit  length  of  cable.  These  two  models 
lead  to  the  following  predictions  for  P.. /P-.  The  relations  given  in 
equations  (33)  and  (34)  were  used  to  eliminate  frequency  from  this 
expression. 


Constant  damping  coefficient 


Constant  damping  ratio 


5 


Jl  Jj. 


Si  f  1 
32  J2 


4V 

L2  U22 


(40) 


P1  L1  U1 

—  ■  -L  (41) 

P2  L2U2 


These  results  suggest  that  in  general  the  response  spectrum  of 
field  cables  exposed  to  shear  flows  will  be  biased  toward  those  modes 
excited  by  regions  of  higher  flow  velocity.  To  be  more  conclusive 
would  require  a  field  experiment  with  a  known  spatially  varying  flow. 
Furthermore,  our  present  knowledge  of  cable  damping  mechanisms  prevents 
more  definitive  conclusions. 


A  Non-Uniform  Cable  in  a  Uniform  Flow 


Consider  a  two  segment  cable  with  lengths,  diameters,  and  virtual 
masses  per  unit  length  L^,  D^,  m^  and  Ls,  Ds>  mg  where  the  subscripts 

l  and  s  refer  to  segments  with  the  larger  and  smaller  diameters,  respec¬ 
tively.  A  uniform  mean  square  lift  coefficient,  C,  is  assumed  to  be 
set  to  an  appropriate  value,  and  Pfc  and  St  are  assumed  to  be  the  same 

for  both  segments.  The  quantity  to  be  predicted  is  Pg/Ps,  the  relative 

spatially  averaged  mean  square  displacement  response,  where  P^  ls  the 
spatially  averaged  mean  square  displacement  response  of  the  modes  ex¬ 
cited  by  flow  incident  on  and  P_  is  the  response  to  flow  incident  on 
D  .  If  the  longest  excited  spatial  wavelength  is  assumed  to  be  much 
shorter  then  both  and  Ls,  then  by  piecewise  integration  techniques 
similar  to  that  in  the  previous  example  the  relative  mean  square  re¬ 
sponse  weighting  P-/P-  may  be  obtained,  as  shown  in  equations  (42)  and 
(43).  1 

Regardless  of  the  chosen  form  of  damping  the  response  is  strongly 
weighted  in  favor  of  the  large  diameter  segment.  The  major  and  unjus¬ 
tified  assumption  in  the  above  example  is  that  there  is  a  uniform  mean 
square  lift  coefficient,  C,  that  ls  the  same  over  both  cable  segments. 

At  the  present  time  we  do  not  have  a  sufficient  understanding  of  the 
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fluid  cable  Interaction  to  asalgn  appropriate  values  for  C  on  each  seg¬ 
ment  of  a  cable  with  varying  diameter.  Nonetheless,  the  result  of  the 
above  analysis  suggests  that  the  cable  response  will  be  highly  dominated 
by  the  modes  excited  by  the  shedding  from  the  larger  diameter  segments 
of  a  cable.  This  is  an  important  qualitative  result  in  tinderstanding 
and  predicting  cable  strumming. 


Constant  damping  coefficient 


Constant  damping  ratio 


- 


fj  S 

Js  Js 


-  5 


is 


_  h  da  !» 

P8  Ls  D4 
8 


(42) 


Lf  Di  "s 

ps  "  Ls  D5  \ 
s 


(43) 


To  test  this  qualitative  prediction  a  simple  experiment  was  con¬ 
ducted  by  slowly  towing  a  nearly  vertical,  heavily  weighted,  composite 
cable  over  the  stern  of  a  ship.  Due  to  length  constraints  on  this  paper 
only  a  brief  account  will  be  presented  here.  A  more  complete  descrip¬ 
tion  may  be  found  in  reference  [2]. 


COMPOSITE  CABLE  EXPERIMENT 


A  two  diameter  composite  cable  was  fabricated  from  a  single  39 
meter  length  of  the  .0071  m  diameter  wire  rope  tested  at  Castine.  The 
rope  had  three  lays  with  an  exterior  sheath  of  polypropylene.  The 
sheath  and  two  of  the  lays  were  removed  from  half  of  the  cable  length. 
This  resulted  in  a  diameter  ratio  D^/D  -  3.0,  where  the  subscripts  1 
and  s  denoted  large  and  small.  A  large  weight  was  attached  to  the  end 
with  the  larger  diameter  D  .  The  cable  was  then  lowered  to  its  full 
length  over  the  stern.  As  the  vessel  moved  slowly  forward  at  speeds 
less  than  0.5  m/s,  the  cross  flow  and  in  line  acceleration  response  at 
several  points  near  the  upper  end  was  recorded.  After  recording  suffi¬ 
cient  data  the  cable  was  reeled  in  3.05  m,  resulting  in  an  Increase  in 
the  ratio  of  immersed  lengths,  L^/L  from  1.16  to  1.49.  Again  the  re¬ 
sponse  of  the  small  diameter  section  was  recorded,  followed  by  pulling 
in  another  3.05  m  section.  This  process  was  repeated  until  only  the 
large  diameter  cable  remained  in  the  water  and  L^/L  •  «.  At  this 
point  the  cable  was  reversed  and  the  experiment  wassrepeated.  This 
time  the  response  on  the  larger  cable  was  measured  for  discrete  values 
of  L^/L  from  .86  to  0.0.  After  computing  corrections  for  mode  shape 
and  mass  loading  of  the  accelerometers  the  average  mean  square  displace¬ 
ment  response  ratios  P^/P  were  computed  and  compared  to  predicted  val¬ 
ues.  The  measured  values  of  p£  and  P  reflect  the  sum  of  both  in  line 
and  cross  flow  spectral  components.  However,  the  cross  flow  response 
was  by  far  the  dominant  component.  The  response  typically  included 
approximately  ten  modes.  Table  III  presents  predicted  and  measured 
values  for  the  experiment  with  the  large  cable  on  top.  The  reverse 
experiment  with  the  small  cable  at  the  top  confirmed  these  results. 

Regardless  of  the  damping  model  chosen  the  results  are  remarkably 
close.  Even  for  the  composite  cable  that  had  the  larger  diameter  cable 
spanning  only  12Z  of  its  total  submerged  length,  over  90%  of  the  mean 
square  displacement  response  was  generated  by  flow  over  the  larger  dia- 


meter  cable  segment.  The  predictions  for  this  example  coae  froa  a 
model  which  was  dalaed  to  be  at  best  a  qualitative  Indicator  of  expec¬ 
ted  response.  The  authors  still  hold  this  to  be  true.  The  predictions 
would  likely  diverge  substantially  from  measured  values  at  very  low 
values  of  L^/CL^  +  Lg). 

The  response  of  a  composite  cable  when  dominated  by  nodes  excited 
in  the  larger  diameter  segment  has  several  important  but  not  necessar¬ 
ily  obvious  features.  First,  each  excited  mode  has  a  mode  shape  which 
spans  the  entire  cable  length.  However,  because  of  the  change  in  the 
mass  per  unit  length,  the  wave  length  and  mode  shape  amplitude  are  lar¬ 
ger  in  the  smaller  diameter  segnent.  In  the  experiment  just  described 
whenever  the  larger  diameter  cable  spanned  more  than  a  third  of  the 
total  length,  the  modes  excited  by  flow  over  the  larger  diameter  con¬ 
tributed  greater  than  99Z  of  the  mean  square  response,  independent  of 
whether  the  response  was  measured  on  the  large  or  small  cable.  When 
measured  on  the  large  diameter  cable  the  response  was  limited  to  approx¬ 
imately  1/3  of  its  own  diameter,  rms.  When  measured  on  the  small  cable 
the  response  typically  exceeded  2  diameters  (of  the  small  cable)  rms. 
This  was  a  necessary  consequence  of  the  relative  amplitudes  of  the  mode 
shapes  in  the  two  segments.  From  laboratory  experiments  with  driven 
cylinders  the  lift  coefficient  has  been  observed  to  decrease  dramatic¬ 
ally  when  the  amplitude  approaches  one  diameter.  This  suggests  that 
under  such  circumstances  the  mean  square  lift  coefficient  on  the  smaller 
segment  drops  to  near  zero,  which  would  explain  the  greater  than  pre¬ 
dicted  dominance  of  the  excitation  on  the  larger  diameter  segment. 

This  also  suggests  that  the  small  diameter  segment  acts  aB  a  damper  or 
a  sink  of  vibrational  energy,  thus  Increasing  the  modal  damping  of  the 
excited  modes. 

TABLE  III.  Prediction  vs  measured  response  for  a  composite  cable. 


Li  +  Ls 

.12 

.24 

.33 

.40 


Predicted  for 


Measured 


'  Vj. 

pi 

p* 

Pi 

Pi  +  Ps 

Pi  +  Ps 

Pi  +  P, 

.71 

.88 

.92 

00 

.94 

.98 

.89 

.96 

.99 

CN 

o\ 

• 

.97 

.999 

CONCLUSIONS 

Cable  strumming  response  for  real  world  cable  systems  in  the  oceans 
is  not  adequately  described  by  resonant  lockin  models.  The  phenomenon 
is  extremely  complex  and  yet  well  described  using  a  stochastic  frame¬ 
work.  The  key  response  characteristics  are: 
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1)  The  mean  square  displacement  response  (normalised  by  cable 
diameter)  is  limited.  This  limit  is  relatively  Independent 
of  cable  length,  virtual  mass,  damping,  number  of  excited 
modes,  and  flow  speed  for  subcritical  Reynolds  numbers. 

2)  The  excitation,  arising  from  vortex  shedding,  is  a  narrow 
band  random  process  centered  about  a  local  "Strouhal"  shed¬ 
ding  frequency. 

A  simple  model  incorporating  the  above  characteristics  explains 
much  of  what  is  experimentally  observed  in  large  ocean  cable  systems. 
The  effect  of  a  spatially  non-uniform  current  was  predicted  and  remains 
to  be  verified  by  experiment.  The  effect  of  a  non-uniform  cable  dia¬ 
meter  was  examined  and  quail a lively  borne  out  by  experiment . 

In  summary  this  initial  attempt  at  a  stochastic  model  has  provided 
considerable  Insight  to  vortex  excited  cable  vibration.  The  model  may 
be  improved.  Additional  work  on  the  mean  square  lift  coefficient  model 
is  warranted  as  well  as  further  research  on  modal  damping  and  its  role 
in  vortex  excited  vibration  response. 
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NOMENCLATURE 

A2: 

B(f,x): 

B1(f): 

B2(f): 

C: 

Cj^x.t): 

C  : 

P 

D(x) : 


mean  square  amplitude  limit  in  cable  diameters 
local  excitation  bandwidth  function 
excitation  bandwidth  function  for  cable  segment  1 
excitation  bandwidth  function  for  cable  segment  2 
spatially  uniform  mean  square  lift  coefficient 
local  instantaneous  mean  square  lift  coefficient 
speed  of  wave  propagation  in  a  cable 
local  cable  diameter 


f(x.t): 


f  : 
P 


f§(x): 


diameter  of  large  diameter  cable  segment  ' 
diameter  of  small  diameter  cable  segment 
local  Instantaneous  force  per  unit  length 
natural  frequency  of  mode  j 
frequency  of  highest  spectral  peak 
local  shedding  frequency 
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H(x) : 

Hj (f) : 

L: 


local  height  of  spectral  box 
transfer  function  for  node  j 
cable  length 

length  of  cable  segment  1 

length  of  cable  segment  2 

length  of  large  diameter  cable  segment 

length  of  small  diameter  cable  segment 

virtual  mass  of  cable;  also  modal  mass 

virtual  mass  per  unit  length  of  large  diameter  cable 

virtual  mass  per  unit  length  of  small  diameter  cable 

spatial  average  mean  square  displacement  response  in 
bandwidth  1 

spatial  average  mean  square  displacement  response  in 
bandwidth  2 


V 

R.(x,y,T): 

Sc(x,y,f): 

Sc(x,f): 

SD(x,y,f): 

Sf(x,y,f): 

Sp  (f): 
PJk 

sv  <f>: 

c  j 

V 

U(x): 


Jr 

h: 


r' 

W(x): 
Af 


5dB' 


average  mean  square  displacement  response  due  to  large 
diameter  cable  excitation 

average  mean  square  displacement  response  due  to  small 
diameter  cable  excitation. 

excitation  bandwidth  fraction 

space-time  correlation  function  for  the  lift  coefficient 
space-time  spectrum  for  the  lift  coefficient 
space-time  spectrum  approximation  for  the  lift  coefficient 
space-time  spectrum  for  cable  displacement  response 
space-time  spectrum  for  excitation  force  per  unit  length 
modal  force  cross-spectrum  ( j ,k) 

modal  force  auto-spectrum  approximation 
"Strouhal"  number  for  a  "non-lockin"  cable 
local  normally  incident  flow  velocity 
normally  incident  flow  velocity  in  region  1 
normally  incident  flow  velocity  in  region  2 
reduced  velocity 
local  width  of  spectral  box 
bandwidth  of  5dB  spectral  band 


6(x-y):  derac  delta  or  impulse  function 


Kronecher  delta 

:  modal  damping  ratio  for  mode  j 
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p  :  density  of  the  fluid 

^ (x) :  eigenfunction  for  node  j 
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